Vacuum polarization of the quantized massive scalar field in Reissner-Nordstrom 

spacetime 

Jerzy Matyjasek, Dariusz Tryniecki and Katarzyna Zwierzchowska 

Institute of Physics, Maria Curie- Sklodow ska University 

pi. Marii Curie- Sklodowskiej 1, 20-031 Lublin, Poland 

(Dated: May 11, 2010) 

The approximation of the renomialized stress-energy tensor of the quantized massive scalar field 
in Reissner-Nordstrom spacetime is constructed. It is achieved by functional differentiation of the 
first two nonvanishing terms of the Schwinger-DeWitt expansion involving the coincidence limit of 
the Hadamard-Minakshisundaram-DeWitt-Seely coefficients [as] and [04] with respect to the metric 
tensor. It is shown, by comparison with the existing numerical results, that inclusion of the second- 
order term leads to substantial improvement of the approximation of the exact stress-energy tensor. 
The approximation to the field fiuctuation, {(f> ), is constructed and briefiy discussed 
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C^ ■ I. INTRODUCTION 
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Recently, there has been a renewal of interest in calculations of the coincidence Umits of the covariant derivatives 

of various bitensors and bispinors, such as the world function a(x,x'), van Vleck- Morette determinant, bivectors 

I——,' and bispinors of the parallel displacement and the objects constructed form them, with the special emphasis put 

O , on the diagonal Hadamard-Minakshisundaram-DeWitt-Seely (HMDS) coefhcients [an{x,x')] and the renormalized 

^ '■ stress-energy tensor ^^ . This has largely been stimulated by new methods of computations as well as improvement 

?H i of the computer algebra and the recent findings substantially extended previous results p-[l5| • 

^*J- In this article we construct the approximation to the renormalized stress- energy tensor of the massive scalar 

I field with arbitrary curvature coupling in the spacetime of the Reissner-Nordstrom black hole. This tensor has been 

T— I . the subject of previous studies both in the massless and massive cases. The extensive numerical calculations have 

^ ■ been reported in Ref. [16J, where, additionally, the first-order approximation in the large mass limit (which consists 

C , of the terms proportional to m~^) has been developed. These calculations were based on the sixth-order WKB 

^^ ' approximation of the solutions of the scalar field equations and the summation of the thus obtained mode functions. 

^" ] The latter result has been reconstructed with the aid of the Schwinger-DeWitt approximation of the renormalized 

• ■ effective action, Wr, and subsequently generalized to the spinor and vector fields in the R = geometries in Ref. [l|. 

7%^ [ The approximate stress- energy tensor in a general background geometry has been constructed in Ref. [2] • A careful 

pZ. • analysis carried out in Ref. 17 1 demonstrated that the approximation is reasonable for Mm > 2 and should become 

I ' increasingly accurate as the ratio Xc/C decreases, leaving, however, room for further improvement. Here M is the 

• • . mass of the Reissner-Nordstrom black hole, Ac is the Compton length associated with the massive field and C is the 

. ^ ' characteristic length of the background geometry. 

k> [ Our aim is to provide a better approximation than those proposed in Refs. [l|, [3, [3 and this is achieved by the 

5^ ■ inclusion of all relevant terms of the background dimensionality 8 (which equals twice the order of the coefficient [a„] or 

Cu [ the total number of derivatives of the metric tensor in each term) in the renormalized effective action. Such terms are 

proportional to m~'^ and constitute the next-to-leading order of the approximation. The approximate stress-energy 

tensor can be obtained by the functional differentiation of the thus constructed effective action with respect to the 

metric tensor. 

The basic building blocks of the renormalized one-loop Wr are the coincidence limits of the HMDS coefficients 
which are local quantities constructed from the Riemann tensor, its covariant derivatives and contractions. The spin 
of the field and the type of the curvature coupling enters Wr through the numerical coefficients and the background 
dimensionality of [03] and [a^] is 6 and 8, respectively. In general, the coefficient [a„] is a linear combination of the 
Riemann monomials and belongs to ®" j^ T^2n qi where TV^ „ is a vector space of Riemannian polynomials of rank r 
(the number of free tensor indices), degree q (number of factors) and order s (number of derivatives). 

The calculation of the functional derivatives of the effective action with respect to gab is rather tedious and time 
consuming process. Fortunately, for the spherically- symmetric geometries one can considerably simplify calculations. 
Indeed, substituting the line element of the general static and spherically-symmetric spacetime expressed in the 
Schwarzschild gauge into the effective action and performing simple integrations one obtains a reduced functional that 
depends on the two metric potentials. Two components of the stress-energy tensor are given by appropriate Lagrange 
derivatives of the Lagrangian of the reduced action functional with respect to the time and radial components of the 
metric tensor. This approach is justified by the symmetric criticality theorems of Palais [T^ [l^l and the remaining 



components can easily be calculated form the covariant conservation equation. 

The coefficients [a„] are also the basic building blocks of the approximate field fluctuation (0^) and the knowledge of 
[^2], [03] and [04] allows for detailed analysis of the role played by the next-to- leading and the next-to-next-to- leading 
terms. One expects that some general features exhibited by {(fp') are also shared by the stress-energy tensor and if so 
it would be a fortunate circumstance. 

The paper is organized as follows. In Sec. II the basic building blocks of the approximation, i.e., the coincidence 
limits of the HMDS coefficients [03] and [04] as calculated within the framework of the covariant DeWitt method 
are presented in maximally condensed form. The renormalized expectation value {(jj^) in the Reissner- Nordstrom 
geometry is computed in Sec. III. The effective action and the stress-energy tensor of the quantized massive scalar 
fields in the Reissner-Nordstrom geometry is constructed and discussed in Sec. IV. A comparison with the numeric 
results indicate that inclusion of the next-to- leading term substantially improves approximation. 

II. HADAMARD-MINAKSHISUNDARAM-DEWITT-SEELY COEFFICIENTS [03] AND [a^]. 



It is a well-known fact that for sufficiently massive quantized fields, i.e., when the Compton length Ac is smaller 
than the characteristic radius of curvature of the background geometry, the asymptotic expansion of the effective 
action in powers of m~'^ may be used to describe various physical phenomena. It is because the nonlocal contribution 
to the effective action can be neglected whereas the vacuum polarization part is determined by the local geometry. 
The renormalized effective action constructed within the framework of the Schwinger-DeWitt approximation for the 
quantized massive scalar field satisfying the covariant Klein-Gordon equation with the curvature coupling ^ 
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can be written in the form 
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where [a„] is constructed from the Riemann tensor, its covariant derivatives up to 2n — 2 order and app ropriate 
contractions. For the technical details of this approach the reader is referred, for example, to Refs. [20, l2l[ and the 
references cited therein. 

Inspection of Eq. ([2]) shows that the lowest term of the approximate Wr is to be constructed from the (integrated) 
coincidence limit of the fourth HMDS coefficient, [03], whereas the next to leading term is constructed form [04]. Here 
we will confine ourselves to the first two terms of the expansion ([2]) . 

The coefficients an{x,x') satisfy the equation 
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with the boundary condition ao(a;, a;') = 1. Here /S.{x,x') is the van Vleck-Morette determinant and the biscalar 
cr(a;,a;') is equal to one half the square of the distance along the geodesic between x and x' . From the recursive 
relation ^ it is clear that to construct, say, [04], one has to calculate coincidence limits of 03, a^.^i^ and a-^.^i^ir^, which, 
in turn, require calculation of [02] to [a2;ii...i4], and so forth. Both [03] and [04] used in this paper have been calculated 
within the framework of the manifestly covariant method proposed by DeWitt with the aid of the FORM and its 



multithread version TFORM [22|, |23 
with the aid of the package INVAR 



Further simplifications, after appropriate syntax conversion, were carried out 

Since the coefficients [03] and [04] are the basis of the present calculations we shall display them at length. The 
coefficient [03] , when simplified with the aid of the INVAR, can be written in the form 
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Before we proceed to the coefficient [04] iet us discuss brieffy this resuit. The package INVAR tries to expand each 
Riemann monomiai in the basis of the independent Riemann invariants with no free indices {r — 0). This is achieved 
by defining polynomial relations between dependent monomials and the basis. Additionally, a number of dimensionally 
dependent identities have been implemented. All this is of great practical importance since it is common that small 
changes in the computational strategy can yield great differences in the theoretically equivalent results and the explicit 
demonstration of their equality is extremely difficult. Having at one's disposal a basis and rules provided by INVAR 
one can easily establish equivalence of the results by construction a unique set of numerical coefhcients and our 
calculations of [03] and [04] have been compared and checked this way. It should be noted that the number of terms 
in [as] has been reduced as compared to the result presented, e.g., in Ref. '26i|, but, of course, they are identical up 
to relatively simple identities satisfied by the Riemann tensor. 

The coefficient [04] is, on the other hand, extremely complicated and even after massive simplifications it consists 
of 113 terms of dimensionality of [1/length]^. The coefficient [04] can be written in the form 

[a,] = af + e4'^ + eaf^ + C^af + ^af^ , (9) 
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Although the total divergences can be discarded when working with the effective action we have retained all the terms 
both in [03] and [04], simply because they are interesting in their own right and can be used in further calculations 
of, for example, the field fluctuation, {(f)^). 

III. FIELD FLUCTUATION 

Before proceeding to the calculation of the renormalized stress-energy tensor in the Reissner-Nordstrom geometry 
let us construct {(fP') as it shares some of the general features of the full stress-energy tensor while simultaneously 
being calculationally less involved. For example, if the next-to-leading term in the expansion leads to substantial 
improvement of the result it is likely that the quality of the analytic approximation of the stress-energy tensor would 
also improve. Similarly, if the vacuum polarization diverges on the event horizon it is quite probable that the stress- 
energy tensor is also divergent there. 

The vacuum polarization has been studied by a number of authors both in the massive and massless cases. (See, 
e-g. [23, [23 and the references cited therein). Generally, within the Schwinger-DeWitt framework the field fiuctuation 
can be expressed in terms of the coincidence limit of the coefficients [a„] and is given by 
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FIG. 1: This graph shows the rescaled values of the vacuum polarization (0^)1 [A = 16tv^AI^] at the event horizon for the 



massive scalar field with mM ■ 
corresponds to {(f>'^)2- 



2 The solid line corresponds to (</> )4 the dashed line corresponds to 



and the dotted line 



where fc — 1 is the number of terms retained in the expansion. The leading term of the expansion is to be constructed 
from the coefficient [02] , which, for the scalar field satisfying Eq. ([T]) is given by 
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Having at our disposal compact expressions describing the first three coefficients of the expansion (|15p we can analyze 
the influence of the higher order terms on the final result. Routine calculations carried out in the Reissner- Nordstrom 
geometry give 
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Now, in order to gain insight into the nature of the thus constructed approximations let us compare {41'^) 2 (the leading 
term), {(j)'^)^ (((/>^)2 plus the next-to-leading term) and {(j)'^)i ((0^)3 plus the next-to-next-to-leading term). In Figs. 



1 and 2 



for fc = 2,3,4 at the event horizon is plotted against the admissible values of \Q\/M. It is seen that 



the second-order term considerably modifies the main approximation. On the other hand, however, ((/>^)3 and {(t>''^)4 
differ only slightly and consequently the changes in the field fluctuation caused by [04] are relatively small. Similarly, 
for a given Q the modiflcation caused by the next-to-next-to-leading term is small everywhere outside the event 
horizon. Therefore, it is reasonable to retain only the first two terms of the expansion (|15p. In what follows we 
shall demonstrate that a similar pattern holds for the stress-energy tensor and the first two terms of the expansion 
([2]) provide a good approximation. There is, however, a profound difference between the two objects: to evaluate 
the renormalized stress-energy tensor on has to retain the terms constructed from [03] and [04] whereas analogous 
calculations of the field fiuctuation require [02] and [a^]. 
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FIG. 2: This graph shows the rescaled values of the vacuum polarization {<^)i [A = IGtt^M^] at the event horizon for the 
massive scalar field with mM = 2. The solid line corresponds to {<?!>^)4 the dashed line corresponds to ( 



and the dotted line 
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IV. THE APPROXIMATE STRESS-ENERGY TENSOR 



A. Effective action 



Having at one's disposal the approximate effective action, W^j, the stress-energy tensor can be obtained form the 
standard formula 



y/gogab 

The total action that leads to the scmiclassical Einstein field equations can be written in the form 



Stotai = Y^ / Rg'^^d'^x + 5™ + Wr, 



where Sm is the action of the classical sources and 
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Now, in order to calculate the functional derivatives of the action with respect to the metric tensor the maximally 
simplified coefficients [03(0;, x')] and [a4{x,x')] have once again been converted into the FORM language. This allows 
us to make use of our powerful FORM codes for functional derivatives of the effective action. The final (mildly 
simplified) result for the renormalized stress-energy tensor of the massive scalar field with an arbitrary curvature 
coupling consists of a few thousand rather complicated terms and is stored in the FORM format. It will not be 
displayed here because it is too bulky. All our calculations in various geometries are always carried out using this very 
tensor and the results are compared (if necessary) with the analogous results constructed using different (possibly 
simpler) computational strategies. 

For the spherically symmetric line element expressed in the Schwarzschild gauge 
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one can save a lot of work by using the reduced action functionals. Indeed, since the coefficients [a„] are constructed 
from various curvature invariants they depend solely on the functions /^"'(t") and h^^'{r), their radial derivatives and 
the radial coordinate. Therefore, one can easily perform simple integrations and reduce the problem to variations 
with respect to the functions /'•^•'(r) and h'-^^r). 

The reduced action functional of the quantum part of the total action is 
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Although stih tedious, this method requires substantiahy smaher number of operations than the general one. This 
procedure yields T" and T'''' components of the renormalized stress-energy tensor; the third algebraically independent 
component, Tg — Tf, can be obtained from the covariant conservation equation. 



The quantum part of the total Lagrangian can schematically be written as 
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/('^^ and /i'*''-' denote a fe— th derivative of f'^^^r) and ft,'"^(r), respectively. Note that the coefficients {32TT^m^^)~-^ 
have been absorbed into the definition of Ci . 

The full form of £1 and £2 calculated for the line element (IMl) . when expanded, consists of 531 and 2157 terms, 
respectively, and will not be presented here for obvious reasons. In practice, there is no need to retain all the terms 
in Eqs. Q and ©. For example, for the R ^ class of metrics which is our main interest here all the terms in Eqs. 
([7]) , (HI) and (|12m4p do not contribute to the final result and the number of terms in 
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as well as ai and 
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is substantially reduced. Further simplifications can be obtained by neglecting the total divergences. 



B. Approximate stress-energy tensor in Reissner-Nordstrom geometry 

Now the stress-energy tensor can be obtained from the Euler-Lagrange equations 
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where n(i) and s(i) can easily be inferred form Eqs. ([25)1 and (P7)) . The angular components can be easily obtained 
from the covariant conservation equation V aT°'^ — 0, which, for the line element (l24l) . reduces to 
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where i = 1, 2. 

Before we start calculations in the Reissner-Nordstrom geometry let us discuss some features of the approximation 
to the stress-energy tensor which can be deduced from the coefficients [03] and [04]. First observe that the stress-energy 
tensor depends linearly on the coupling parameter S, and can be written as a sum of tensors of the type 
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where r] = ^ — 1/6. Further, it should be noted that since T^ is constructed solely from the Riemann tensor and its 
covariant derivatives it is regular for regular metrics. This property is guaranteed by the polynomial character of the 
result and the factorization: 
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where P^'^'{r) are the regular functions. 
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Although the first-order approximation to the stress-energy tensor in the Reissner-Nordstrom geometry is known, 
we shall display it for reader's convenience. Making use of our general formulas, after some algebra, we conclude that 
the approximation has the form p2|) . where 
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This tensor is identical to that constructed in Refs. |il.ll6|. 

Now, let us return to the tensor Ta ■ The calculations that lead to this object are far more complicated than the 
analogous calculations of the first-order term and require heavy use of the computer algebra. All these efforts however 
will pay off and give us substantially better approximation. In recent publications [J, |29[ it has been argued that 
the minimal approximation constructed within the framework of the Schwinger-DeWitt method should consist of the 
two first terms of the expansion ([2]). This observation was based on the analyses carried out in the Schwarzschild 
and the Bertotti-Robinson geometries. Specifically, it has been demonstrated that the approximation of the stress- 
energy tensor in the Schwarzschild spacetime constructed form [as] and [04] is substantially better that the analogous 
approximation calculated from the coefficient [03] alone and this by itself justifies the introduction of the second order 
term in that case. We shall show that similar behavior occurs in the spacetime of the Reissner- Nordstrom black hole. 
Moreover, the higher order terms may dramatically change the type of the solutions of the semiclassical Einstein field 
equations. An interesting example in this regard is given by the Bertotti- Robinson geometry [30, |3ll- Specifically, it 
can be shown that although the Bertotti-Robinson geometry is a self-consistent solution of the semiclassical Einstein 
field equations with the source term given solely by the leading term of the renormalized stress-energy tensor [1|, i32h34j 
it does not remain so when the next-to-leading term is taken into account. To guarantee that the Bertotti- Robinson 
spacetime is the solution of the semiclassical equations one has to introduce the (negative) cosmological constant. It 
should be noted that addition of the electric charge to the system does not change this behavior. 

Thee second-order term has, as expected, the form ([5^. where 
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FIG. 3: This graph shows the rescaled values of Co [A = 90(8M) tt ] as function of a: = (r — r+)/M for massive scalar field 
with mM = 2 and \Q\/M — 0.95. The solid line corresponds to the improved approximation whereas the dashed line to the 
first order Schwinger-DeWitt approximation. In the small panel the near horizon behavior of Cg is displayed 
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The angular components of the renormalized stress-energy tensor of the massive scalar field in the Reissner- Nordstrom 
geometry has been calculated numerically for mM = 2 and \Q\/M = 0.95. (Similar calculations have also been 
carried out for \Q\ — 0.) However, in view of further applications, it is preferable to have at one's disposal simple and 
accurate general analytic formulas describing the functional dependence of the stress-energy tensor on the metric. Such 
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FIG. 4: This graph shows the rescaled values of Dg [A = 90(8M)*7r'^] as function of x = (r — r+)/M for massive scalar field 
with mM = 2 and \Q\/M — 0.95. The solid line corresponds to the improved approximation whereas the dashed line to the 
first order Schwinger-DeWitt approximation. In the small panel the near horizon behavior of Dg is displayed. 

general formulas can easily be applied to the concrete line element provided some general requirements concerning the 
geometry and the mass of the quantized field are satisfied. Consequently, it is of interest to compare the approximation 
constructed from [03] and [04] with the results of the numerical calculations of the conformal and nonconformal 
contribution to the total stress-energy tensor 
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as presented in ref. [lal 

A comparison of Figs. 3 and 4 of the present paper with the Figs. 10 and 11 of Ref. [16|] clearly shows that 
although the first order approximation correctly reproduces qualitative behavior of C| and Dg , the inclusion of the 
next-to-leading term substantially improves the approximation of the stress-energy tensor even in the closest vicinity 
of the event horizon. One expects that this approximation is even better for mM > 2. A lesson that follows from this 
demonstration is that the next-to-leading term plays, or at least may play, an important role in the calculations and 
it can be ignored only after careful examination. Similar behavior of the next-to-leading term of the stress-energy 
tensor in the Schwarzschild spacetime has been reported in Ref. t4| . One expects therefore that this pattern holds for 
all values satisfying < \Q\/M < 0.95. There is also good reason to believe that it is true for all admissible values of 

q- 

At the event horizon of the extreme Reissner-Nordstrom black hole the stress-energy tensor can be written in the 
remarkably simple form 
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As the geometry in the vicinity of the degenerate Reissner-Nordstrom solution is precisely that of the Bertotti- 
Robinson one can easily calculate the approximate stress-energy tensor. Indeed, due to homogeneity {Rabcd;e = 0), 
vanishing of the Ricci scalar and the Weyl tensor {Cabcd = 0), the stress-energy tensor can be expressed solely in 
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terms of the Ricci tensor. Making use of the general formulas the first order approximation to the stress-energy tensor 
can be written in the form 

1260 «" ' V ; 

whereas when the above conditions are satisfied the next-lo-leading term reads 

75600 ''b cd ij \ J 

Simple calculation shows that (H51) and (|i^ in the Bertotti-Robinson geometry are precisely equivalent to the stress- 
energy tensor at the degenerate horizon of the extreme Reissner-Nordstrom black hole and this may be regarded as 
the additional useful check of the calculations. 

V. FINAL REMARKS 

In this work our goal was to construct the approximate field fluctuation and renormalized stress-energy tensor of the 
quantized massive field in the spacetime of the Reissner-Nordstrom black hole and to investigate how the higher-order 
terms of the expansions ^ and p^ affect the final results. The general formulas describing the both quantities 
are extremely complex, but, fortunately, there are massive simplifications when applied to the static and spherically 
symmetric geometries. A comparison with the numeric calculations reported in a classic paper by Anderson, Hiscock 
and Samuel |16| shows that the next-to-leading term substantially improves the approximation. It has been found 
that in both cases the minimal approximations are to be constructed from the first two terms of psp and ^ for 
the field fluctuation and the stress-energy tensor, respectively. Although we have constructed the general form of the 
stress-energy tensor up to the next-to-leading terms by functional differentiation of the action functional with respect 
to the metric, here we proposed a computationally simpler method in which the reduced action functionals are varied 
with respect to the functions guii") and grrir). Both methods, when overlap, give, of course, identical results. 

We hope that our results will be of use in further calculations. We indicate a few possible directions of investigations. 
First, it would be interesting to analyze the back reaction of the quantized massive field upon the geometry of the 
Reissner-Nordstrom black hole. Due to simplicity of the stress-energy tensor in the Reisssner-Nordstrom spacetime 
the quantum-corrected metric can easily be constructed performing two elementary quadratures. Since the quantum 
part in the right hand side of the semiclassical Einstein field equations is calculated in a large mass limit it is 
purely geometric quantity and can be expressed solely in terms of the Riemann tensor, its covariant derivatives and 
contractions. This allows to treat the semiclassical theory as the higher derivative theory and construct various 
characteristics encoded in the geometry of the quantum-corrected black hole such as location of the horizons in 
nondegenerate as well as degenerate case [33, llg, equations of motion of the test particles fSg, |33|, temperature 
and entropy |37h42| | and energy- momentum complexes [43| . Further, construction of the approximate stress-energy 
tensor as well as the field fluctuation in more complex backgrounds and the accompanying numerical calculations 
would certainly strengthen our understanding of the problem. Especially interesting in this regard is the problem 
of the lukewarm [28[ and ultraextremal [4J] black holes. Finally, an important and interesting continuation of the 
calculations presented in this paper would be construction of the of the next-to-leading term of the spinor and vector 
fields. We intend to return to this group of problems elsewhere. 



[11 


J. 


[21 


J. 


|3| 


R 


|4| 


J. 


[51 


A 


|6| 


B 


|7| 


T 


|8| 


P. 


[91 


P. 


[101 


L. 


[111 


I. 


|12| 


I. 


[131 


I. 



Matyjasek, Phys. Rev. D61, 124019 (2000). 
Matyjasek, Phys. Rev. D63, 084004 (2001). 

T. Thompson and J. P. S. Lemos Phys. Rev. D80, 064017 (2009). 
Matyjasek and D. Tryniecki, Phys. Rev. D 79, 084017 (2009). 
. C. Ottewill and B. Wardell (2009), arXiv:0906.0005 [gr-qc] . 
S. DeWitt, Dynamical Theory of groups and fields (Gordon and Breach, New York, 1965). 
Sakai, Tohoku Math. J. 23, 589 (1971). 
B. Gilkey, J. Differential Geometry 10, 601 (1975). 
B. Gilkey, Trans. Am. Math. Soc. 225, 341 (1977). 

Parker, in Recent Developments In Gravitation (Academic Press, New York, 1979). 
G. Avramidi, Ph.D. thesis, Moscow State University (1986). 
G. Avramidi, Theor. Math. Phys. 79, 494 (1989). 
G. Avramidi, Nucl. Phys. B355, 712 (1991). 



12 



[14] 


P. 


[15] 


A. 


[16] 


P. 


[17] 


B. 


[18] 


R. 


[19] 


S. 


[20] 


A. 


[21] 


V. 


[22] 


J. 


[23] 


M 


[24] 


J. 


[25] 


J. 


[26] 


L. 


[27] 


V. 


[28] 


E. 


[29] 


J. 


[30] 


I. 


[31] 


B. 


[32] 


L. 


[33] 


V. 


[34] 


J. 


[35] 


A. 


[36] 


M 


[37] 


J. 


[38] 


R. 


[39] 


V. 


[40] 


M 


[41] 


M 


[42] 


J. 


[43] 


J. 


[44] 


J. 



Amsterdamski, A. L. Berkin, and D. J. O'Connor, Class. Quant. Grav. 6, 1981 (1989). 

E. M. van de Ven, Class. Quant. Grav. 15, 2311 (1998). 

R. Anderson, W. A. Hiscock, and D. A. Samuel, Phys. Rev. D 51, 4337 (1995). 

E. Taylor, W. A. Hiscock, and P. R. Anderson, Phys. Rev. D61, 084021 (2000). 

S. Palais, Comm. Math. Phys. 69, 19 (1979). 

Deser and B. Tekin, Class. Quant. Grav. 20, 4877 (2003). 

O. Barvinsky and G. A. Vilkovisky, Phys. Rept. 119, 1 (1985). 

P. Frolov and A. I. Zelnikov, Phys. Rev. D29, 1057 (1984). 
A. M. Vermaseren (2000), math-ph/0010025. 
. Tentyukov and J. A. M. Vermaseren (2007), hep-ph/0702279. 

M. Martin-Garcia, D. YUanes, and R. Portugal, Comp. Phys. Commun. 179, 586 (2008). 
M. Martin-Garcia, R. Portugal, and L. Manssur, Comp. Phys. Commun. 177, 640 (2007). 
Parker and D. J. Toms, Phys. Rev. D31, 953 (1985). 

P. Frolov and I. D. Novikov, Black hole physics (Kluwer Academic Publishers, Dordrecht, 1998). 

Winstanley and P. M. Young, Phys. Rev. D77, 024008 (2008). 
Matyjasek and D. Tryniecki, Mod. Phys. Lett. A24, 2517 (2009). 
Robinson, Bull. Acad. Pol. Sci. 7, 351 (1959). 

Bertotti, Phys. Rev. 116, 1131 (1959). 

A. Kofman and V. Sahni, Phys. Lett. 127B, 127 (1983). 

Sahni and L. A. Kofman, Phys. Lett. A117, 275 (1986). 
Matyjasek and O. B. Zaslavskii, Phys. Rev. D64, 104018 (2001). 

Dobado and A. L. Maroto, Phys. Lett. B316, 250 (1993). 
. Lu and M. B. Wise, Phys. Rev. D47, 3095 (1993). 
Matyjasek, M. Telecka, and D. Tryniecki, Phys. Rev. D73, 124016 (2006). 

M. Wald, Phys. Rev. D48, 3427 (1993). 

Iyer and R. M. Wald, Phys. Rev. D50, 846 (1994). 
. Visser, Phys. Rev. D48, 5697 (1993). 
. Visser, Phys. Rev. D48, 583 (1993). 
Matyjasek, Phys. Rev. D74, 104030 (2006). 
Matyjasek, Mod. Phys. Lett. A23, 591 (2008). 
Matyjasek and O. B. Zaslavskii, Phys. Rev. D71, 087501 (2005). 



